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ABSTRACT 

We present new calculations of transit spectra of super-Earths that allow for atmospheres with 
arbitrary proportions of common molecular species and haze. We test this method with generic 
spectra, reproducing the expected systematics and absorption features, then apply it to the nearby 
super-Earth GJ 1214b, which has produced conflicting observational data, leaving the questions of a 
hydrogen-rich versus hydrogen-poor atmosphere and the water content of the atmosphere ambiguous. 
We present representative transit spectra for a range of classes of atmosphere models for GJ 1214b. Our 
analysis supports a hydrogen-rich atmosphere with a cloud or haze layer, although a hydrogen-poor 
model with ^10% water is not ruled out. Several classes of models are ruled out, however, including 
hydrogen-rich atmospheres with no haze, hydrogen-rich atmospheres with a haze of ~0.01 /xm tholin 
particles, and hydrogen-poor atmospheres with major sources of absorption other than water. We 
propose an observational test to distinguish hydrogen-rich from hydrogen-poor atmospheres. Finally, 
we provide a library of theoretical transit spectra for super-Earths with a broad range of parameters 
to facilitate future comparison with anticipated data. 



1. INTRODUCTION 

Ground- and space-based transit surveys are discover- 
ing a growing number of exoplanets in the super-Earth 
range, defined alternately as between 1 and 10 Earth 
masses or between 1 and 2 Earth radii. The largest such 
survey, being conducted by the Kepler spacecraft, has 
already announced hundreds of super-Earth candidates 
(Kepler Objects of Interest) a nd one con firmed sviper- 
Earth in the "habitable zone" (jBorucki e t al. 2011), al- 
though only a few of these candidates have been con- 
firmed. Nevertheless, these surveys have yielded several 
exoplanets for which it is possible to begin character- 
izing their atmospheres via transit spectroscopy. This 
has been attemp t ed by several groups for GJ 1214b 
(IBean et all [Mil IBerta et al. 2011b; CroU et al. 2011; 
Crossfield. Barman, fc Hansen .2011.: .Desert et al.. .2011.: 
De Mooii et al.l 1201 If) , an estimated 6.55 Earth- mass 
(M0) planet orbiting an M-d warf star with an equilib - 
rium temperature near 500 K (jCharbonneau et al.ll2009t ). 
with conflicting results. In particular, measured transit 
depths in the K and Kg bands are in direct conflict, and 
proposed atmospheric models have been unable to flt ob- 
servations in other bands simultaneously. 

GJ 1214's faintness {mv — 14.67) makes it a challenge 
for spectroscopic observations. However, the recent dis- 
covery of transits of 55 Cue e around a naked-eye star 
(Winn et al. 2011) holds promise for improved spectro- 
scopic characterizations in the near future. These stars 
bring their own challenges of smaller transit depths and 
fewer nearby comparisons stars, but efforts are nonethe- 
less underway to characterize them. 

While transit and radial velocity observations can 
together constrain a planet's mass and radius (and 
thus its density), this leaves signiflcant degeneracies 
in its bulk composition. Planets with similar masses 
and radii can take on different compositions, ranging 
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from a mini-Neptune with a thick H/He layer to a 
true super-Earth composed mostly of silicates with a 
high-molecular-weight atmosphere to a "water world" 
whose bulk compositi on and atmosphere both have 
high water content iV alencia. Sasselo v. fc O'Connell 
2007t iFbrtnev. Marlev. fc Barnes 2007^ iSeager et al 



20071 iSotin. Grasset. fc Mocauetl 120071: iRogers fc Seager 
2010). To further confuse matters, recent observations 
suggest that some giant exoplanets may have signifl- 
cantly higher C/0 ratios than their parent stars, opening 
up the additional pos sibility of carbon-rich compositions 
([Madhusudhan et al.ll20l0| ). 

GJ 1214b faces the additional complication that while 
its measured mass and density suggest a hydrogen- 
dominated atmosphere, the measured radius of the par- 
ent star, GJ 1214, is sig nificantly larger than pre - 
dicted by stellar models (jCharbonneau et al.l l2009| ). 
If the star is 15% smaller than was measured, as 
these stellar models predict, GJ 1214b would be cor- 
respondingly smaller and denser (given the same tran- 
sit depth), and a hydrogen- rich atmosphere would not 
be indicated. This uncertainty means that we cannot 
definitively distinguish between a hydrogen-rich and a 
hydrogen-poor atmosph ere w ithout spectroscop ic data 
(|Miller-Ricci fc Fortnevl[20Tol: iCarter et al.ll201lb . 

Spectroscopic observations of both secondary eclipses 
("emission") and transits ("transmission") of exoplan- 
ets have proven successful in char acterizing the at 
mosp hercs of gas giants , e.g. (iBurrows fc Orto 
20091: iCharbonneauet all 120021: iDeming et al.l 1200 
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Knutson et al.l l2008t iSwain. Vasisht. fc Tinettil l2008l) . 
and observations to characterize super-Earths continue. 
Transit spectra are particularly useful for measuring 
compositions. While secondary eclipses can be used 
to measure temperatures and compositions, the transit 
spectrum is most directly sensitive to the column depth 
at the limb of individual molecular and atomic species 
and is only weakly dependent on temperature (for a given 
scale height). In this paper, we present calculations of 
transit spectra for GJ 1214b for a variety of atmospheres 



in order to improve the current fits to the data. We 
also provide a library of models for future exploration of 
transit spectra of super-Earths. 

Transit spectroscopy measures the apparent size of a 
planet in various wavelengths, which depends on the av- 
eraged properties of the planet's atmosphere. Because 
exoplanets cannot be resolved by current telescope tech- 
nology, this transit spectrum is a measure of the total 
amount of starlight blocked in transit at each wavelength. 
Ignoring the effects of limb darkening, this includes a 
fixed component due to the solid disk of the planet, plus a 
wavelength-dependent component absorbed or scattered 
by the atmosphere. Within the atmosphere, different 
altitudes will contribute different amounts to this com- 
ponent because of the variations with pressure. How- 
ever, molecular absorption features should still emerge 
from these spectra in a predictable way, enabling abun- 
dance characterization of planetary atmospheres based 
on transit observations. Conversely, the absence of these 
features in the observed spectrum might be indicative 
not only of their absence, but of a cloud or h aze layer, as 
is inferred in the atmosphere of HD 189733b ( iPont et aD 
12008). but one should keep an open mind that a small 
scale height that reduces the cross section of the atmo- 
sphere could also cause this effect. The molecular fea- 
tures can be characterized by the relative increase in 
transit depth when compared with the continuum (or 
with a local minimum in transit depth if the continuum 
is obscured), a contrast that we hereafter refer to as the 
"amplitude" of the spectral features. 

If the atmosphere is non-isothermal or not verti- 
cally mixed, this further complicates the results. Also, 
the limb of the planet as seen in transit represents 
the terminator, where there is likely to be circula- 
tion and a change in temperature and composition. 
Indeed, the physics at the terminator is more com- 
plicated than any model to date. However, one- 
dimensional theoretical spectra are not strongly depen- 
dent on these factors and provide valuable zeroth-order 
approximations to estimate atmosph eric compositions 
(jMiller-Ricci. Sasselov. fc Seaged 120091 ) . In our models, 
hydrogen-rich atmospheres are taken to be in chemical 
equilibrium, while hydrogen-poor atmospheres are taken 
to be vertically mixed. 

The physical parameters of an exoplanet should have 
predictable effects on its transit spectrum. These mainly 
act through the scale height of the atmosphere, which is 
given by 

H = ^. (1) 

where k is Boltzmann's constant, T is the atmospheric 
temperature, /z is the mean molecular weight, and g is 
the surface gravity. If the atmosphere is approximately 
isothermal, T is the "equilibrium temperature." 

The amount of starlight absorbed in transit is propor- 
tional to the depth of the atmosphere and, thus, to its 
scale height. This should be visible in the transit spectra 
as proportional changes in the amplitudes of the spectral 
features and the slope of the Rayleigh continuum. There- 
fore, these amplitudes should be directly proportional to 
the temperature of the atmosphere and inversely propor- 
tional to its mean molecular weight. To first order, the 
amplitudes of spectral features should also be inversely 



proportional to surface gravity, but this relation is not 
exact since gravity changes slightly with altitude. In 
very thick atmospheres, the lower gravity at high alti- 
tudes should make the scale height larger, resulting in 
an even larger atmosphere and larger amplitudes than a 
simple estimate might suggest. 

The remaining parameters that affect the transit spec- 
trum are related to the composition structure of the at- 
mosphere. In particular, for a hydrogen-rich atmosphere, 
metallicity is important. While for a hydrogen-rich at- 
mosphere metallicity does not have a strong effect on 
mean molecular weight, the relatively large absorption 
opacities of metal compounds compared with molecular 
hydrogen mean that significantly larger-amplitude spec- 
tral features will be found at larger metallicities than 
at smaller ones. The addition of clouds or haze should 
have the opposite effect. The wavelength dependence of 
extinction opacity for haze is usually smaller than for 
molecular species, so the corresponding spectral features 
should be weaker. In the simplest case of a uniformly 
opaque cloud layer, suggested by CroU et al. (2011) and 
Bean et al. (2011) for GJ 1214b, only molecular absorp- 
tion above the cloud tops contributes to the variation in 
the spectrum, and the molecular features are suppressed. 
All of these relations are seen in our calculated spectra. 

We present new calculations of theoretical transit spec- 
tra of extrasolar super-Earths, including Mie scattering 
and a range of molecular compositions. Modeling of tran- 
sit spectroscopy has been established for super-Earths 
previously, for example, in Miller- Ricci & Fortney (2010), 
and it has been used to try to characterize the atmo- 
sphere of GJ 1214b. However, we present a number of 
new results and models including spectra for sub-solar 
metallicities, enhanced C-to-0 ratios, and explicit calcu- 
lations for a range of haze models. We also outline a 
library of models for super-Earths with different param- 
eters. We present our methodology in Section [51 We 
present general conclusions and demonstrations of the 
systematics described above in Section |3l We apply our 
method to GJ 1214b in Section |4] in an attempt to pro- 
duce a better fit to observations than current models. 
We then present some general conclusions in Section [5l 

2. METHODOLOGY 

For each atmosphere model, we produce a theoreti- 
cal spectrum by computing the transit depth at 2500 
frequency points between 0.3 and 10 /im. At each fre- 
quency point, we numerically integrate the optical depth 
along the line of sight of each impact parameter, then 
integrate the effective area over all impact parameters. 
For transiting planets, we use the initial observed transit 
depth as a fiducial value for the radius and then nor- 
malize the spectrum by computing the average transit 
depth over the observational bandpass, equating this to 
the observed depth. Without independent knowledge of 
the radius of the solid surface, a procedure like this is 
necessary for all models and modelers. 

We use an isothermal temperature-pressure ( T-P) pro- 
file at the equilibrium temperature. However, we set 
T == 470 K for GJ 1214b based on the work of Miller- 
Ricci & Fortney (2010), who find with more ambitious 
calculations that this is approximately correct over the 
pressure range probed by transit spectroscopy, although 
their approach does not fully capture the conditions at 



the terminator and should not be considered more reU- 
able than any other modeL Additionahy, transit spectra 
have been shown not to be strongly dependent on the 
exact T-P profile used fMi ller-Ricci. Sasselov. fc Seaged 
l200 9f). We do account for the variation of gravity with 
altitude when computing the pressure profile, although 
this effect is only a few percent in transit depth, depend- 
ing on the scale height of the atmosphere. We implement 
haze using the Mie approximation for a monodispersed, 
constant abundance haze between upper and lower pres- 
sure levels. We also include the option of completely 
opaque cloud decks by cutting off transmission at all 
frequencies at pressures greater tha n a specified level 
(jCrossfield. Barman, fc Hansenll20TTl ). 

Our chemical equilibrium abundances for molecular 
species were taken from Burrows fc Sharp (1999) and 
molecular line and coUisional opacities were taken from 
Sharp fc Burrows (2007). We compute Rayleigh cross 
sections on a per-molecule basis, based on each species 
electric polarizability and assuming a simple A"'* power- 
law dependence. For haze components, we use Mie scat- 
tering theory for several possible haze species. Poly- 
acetylenes are polymers based on CH units produced by 
UV photolysis of simple hydrocarbons, the simplest be- 
ing (CH)n, and a re suggested to occur on Jupiter an d 
giant exoplanets (jSudarskv. Burrows, fc Hubenvl 120031 ). 
Tholins are more complex heteropolymers produced by 
UV photolysis of hydrocarbons, nitrogen, and oxygen 
and are thought to occur in the atmosphere of Titan 
and on the surfaces o f other outer solar system bod- 
ies (jKhare et al.l |1984[) . They take the form of cross- 
linked chains of a wide variety of organic subunits, in- 
cluding aromatic hydrocarbons, other ri ngs, alkanes, and 
pol ypeptides, with no definite order (jSagan fc Khard 
[l97 9) . We also consider sulfuric acid because of its pres- 
ence on Venus. Complex indices of refraction were taken 
from Khare et al. (1984) for tholins, from Bar-Nun et al. 
(1988) for polyacetylenes, and from Palmer fc Williams 
(1975) for sulfuric acid. Our absorption opacities for pure 
molecular species at a temperature and density represen- 
tative of GJ 1214b are shown in Figure[TJ The extinction 
opacities of hazes of various sizes and compositions are 
shown in Figure [51 

We calculate transit spectra for a variety of atmo- 
spheric compositions. The hydrogen-rich compositions 
we compute include metal abundances varying from 0.01 
to 30 times solar, solar abundances with methane set 
to zero, and otherwise-solar abundances with the abun- 
dances of carbon and oxygen reversed, resulting in a 
C/0 ratio of ^2. We consider an atmosphere with 
no methane because methane is easily destroyed by 
UV photolysis, and it has been suggested as a possi- 
bility to explain the appearance of GJ 1214b's spec- 
trum (iMiller-Ricci fc FoTtaevI [2OIOI : iBean et al.l [20TTI : 
iDesert et al.ll201lD ^ We consider an elevated C/0 ra- 
tio because such a possibility has been advanced for gi- 
ant planets, i ncluding those with a low C/0 ratio in 
the host star (Madhusudhan et al."2010) and for terres- 
trial planets (Kuchncr fc Scagcr 2005). The C/0 ra- 
tio is especially important because it strongly affects 
the abundances of common molecular species such as 
H20,CO,C02,andCH4. 

For hydrogen-poor models, we use vertically mixed at- 
mospheres that can include H2, CO2, H2O, CH4, NH3, 



CO, and N2 in arbitrary proportions by volume. We as- 
sume nitrogen to have Rayleigh scattering opacity only. 

To test the code, we compute transit spectra for atmo- 
spheres with pure molecular compositions and for solar 
metallicities over a range of metallicities, temperatures, 
and opaque clouds placed at varying altitudes. These cal- 
culations successfully reproduce the relations predicted 
by systematics, and they are discussed in detail in Section 
[3l The spectra for pure molecular atmospheres succeed 
in reproducing the molecular opacity features of each 
species (compare Figures [1] and ^ . Smooth regions in 
the spectra represent the Rayleigh continuum (or, if the 
continuum is sufficiently transparent, the cloud tops or 
the planet's solid surface) at wavelengths where absorp- 
tion is negligible. 

We also compute spectra for GJ 1214b with compo- 
sitions similar to those used in Miller-Ricci fc Fortney 
(2010) with a "surface pressure" of 1 bar. At this pres- 
sure level, the atmosphere becomes opaque at all wave- 
lengths in most models, so a deeper atmosphere will 
not affect the spectrum. Our spectra successfully repro- 
duce the molecular features of the Miller-Ricci fc Fortney 
(2010) computations. They have systematically larger 
transit depths than Miller-Ricci fc Fortney (2010) in the 
molecular features by 10%-20%, and they agree to bet- 
ter than 10% with the slope of the Rayleigh continuum. 
The systematic offsets are most likely due to the use of 
different opacity tables. 

3. GENERIC MODEL RESULTS 

For convenience, all generic spectra are computed us- 
ing the physical parameters of GJ 1214b as given in 
Tables [T] and [2l except as specified. They are also 
normalized to the initial measurement of GJ 1214b: 
1.35% (13,500 ppm) in th e MEarth band of 0.7-1.0 ^m 
(|Charbonneau et al.ll2009h . 

Generic spectra computed for hydrogen-rich atmo- 
spheres for which abundances are a fraction or multiple of 
solar are shown in Figure |4l These spectra show transit 
depths increasing with metallicity. This is expected since 
metals are the main source of extinction in a haze-free 
atmosphere. However, for this normalization, the tran- 
sit depth of the Rayleigh scattering tail decreases with 
increasing metallicity. This is because of the normaliza- 
tion to the MEarth band, in which absorption dominates. 
The spectral features in that band increase in amplitude 
with increasing metallicity relative to the Rayleigh con- 
tinuum, so the continuum has a smaller relative transit 
depth at large metallicities. 

Two other hydrogen-rich atmosphere models are shown 
in Figure [5l Swapping the abundances of carbon and 
oxygen from solar levels-in effect, raising the C/0 ratio 
from 0.5 to 2.0-diminishes the spectral features of water 
and enhances those of methane. Some of these features 
overlap and, thus, those wavelengths are not strongly 
affected. The greater extinction opacity of methane at 
near-infrared wavelengths causes the Rayleigh continuum 
to be normalized to a lower transit depth. For the model 
with no methane, the spectrum is more strongly dom- 
inated by water features, and the normalization is not 
strongly affected. 

The effects of temperature on the generic transit spec- 
tra are shown in Figure |6l A higher temperature re- 
sults in a larger scale height, and, as expected, a larger 



continuum slope at short wavelengths. As the tempera- 
ture increases, two major chemical changes are evident 
as well-the suppression of methane and the appearance 
of strong alkali metal lines. 

An opaque cloud layer can suppress the molecular spec- 
tral features if it lies at high enough altitudes. If light is 
not sufficiently blocked as it grazes the cloud tops, the 
observed radius should be very near the cloud tops in 
that wavelength. This effect is demonstrated in Figure 
[71 which includes spectra with cloud layers at different 
pressure levels. Flat, horizontal sections of the spectra 
represent wavelengths where the atmosphere is transpar- 
ent down to the cloud tops. If the clouds are deeper 
than the 1 bar level, they are usually obscured at all 
wavelengths, so higher pressures have no effect on the 
spectrum. However, for some atmospheres composed of 
pure molecular species, such as carbon dioxide, light can 
penetrate deeper at some infrared wavelengths, where 
both absorption and scattering are negligible, so these 
flat regions due to clouds can persist to very high pres- 
sures in such an atmosphere. 

Models for pure atmospheric compositions are shown 
in Figure [31 to be compared with Figure [l] which depicts 
the absorption opacities for the same species at a fixed 
temperature and density. Nitrogen is assumed to have 
Rayleigh scattering opacity only, as are the other species 
at wavelengths outside their respective absorption bands. 
All of these results are consistent with the expected de- 
pendencies on physical and atmospheric parameters and, 
thus, provide a further validation of our methodology. 

We add haze to the atmosphere using the Mie ap- 
proximation for spherical particles. We find that dense 
hazes (n > 10^ cm~^ for most typical sizes) are essen- 
tially opaque in the wavelength range we study and re- 
sult in fiat regions in the spectrum. However, very thin 
hazes (n < 10^ cm~'^ for typical particle sizes) become 
transparent at some wavelengths. In general, these thin 
hazes tend to produce sloped regions of the spectrum 
similar to those produced by Rayleigh scattering. How- 
ever, these regions do not span a wide range of transit 
depths because many molecular features persist to much 
lower pressures than typical hazes. Examples of spectra 
with haze included are shown in Figure [Sj 

4. GJ 1214B 
4.1. Prior Observations 

The photometric and spectroscopic data for GJ 1214b 
available at this writing are shown in Figure [S] Desert et 
al. (2011), working in the Warm Spitzer 3.6 and 4.5 /im 
bands, observed a flat transit spectrum of GJ 1214b and 
concluded that its atmosphere must have a small scale 
height and, thus, a high molecular weight, probably with 
a large proportion of water. Bean et al. (2010), working 
in the / band, came to the same conclusion, although 
they also suggested that the planet could have a layer of 
clouds or hazes above the 300-mbar level. 

This picture changed when Croll et al. (2011) ob- 
served a deep Ks-hand transit, but shallow J and 
CH40n(1.69/im) band transits, suggesting a hydrogen- 
rich, low-molecular-weight atmosphere, observable in the 
Ks band, that included a haze layer to mute features at 
other wavelengths. The anomalous Kg-hand feature was 
attributed to either water or methane. 



However, Crossfield et al. (2011), based on near-IR 
spectroscopy, argue that a hydrogen-rich atmosphere is 
not favored, and that if a large-scale-height atmosphere 
were present, it must be out of chemical equilibrium (e.g., 
depleted in methane) or muted by a haze or cloud layer. 
Bean et al. (2011) revised their earlier assessment with 
data in the R, J, H, and K bands. These data support 
a high-molecular-weight atmosphere, particularly with 
shallow K-band transits. However, their reported results 
in the R-band show some evidence for a Rayleigh scat- 
tering tail at short wavelengths, with a slope consistent 
with a hydrogen-rich atmosphere. 

The most recent observations continue to give conflict- 
ing results. De Mooij et al. (2011), taking measurements 
in a wide range of visible and infrared bands, obtained 
results that reinforce both the conclusion of both appar- 
ent Rayleigh scattering at short wavelengths in Bean et 
al. (2011) and Croll's deep i^^-band measurement, pro- 
viding further evidence for a hydrogen-rich atmosphere. 
However, the small relative size of the observed Ks- 
hand feature, compared with theoretical predictions for a 
hydrogen-rich atmosphere, suggested an improbably low 
metallicity for the planet, despite orbit ing a star with 
normal metallicity (De Mooii et al.ir2011| ) . Finally, Berta 
et al. (2011b) observed what appears to be a feature 
at ^1.5 fim. with the Hubble Space Telescope's WFC3, 
which could correspond to water, the depth of which is 
consistent with a high-molecular-weight atmosphere. 

To explain these conflicting results, one of the most 
widely postulated models is a hydrogen-rich atmosphere 
with a layer of clo uds or haze near or above the 10 
mbar pressure level ([Crofl et al.ll2011l:lBean et al.ll201lD , 
or above the 10 mbar level in Berta et al. (2011b). 
Haze has also been suggested to explain the weak al- 
kali meta l features in tran sits of the giant planet HD 
189733b (|Pont et al.l l2008[ ). The other alternative is 
a high-n iolecular-weight atmosphere with a h igh water 
'^" al.li2011tlBerta et al.ll20liB . 



content ([Desert et 



4.2. Models without Haze 

Our full list of atmosphere models for GJ 1214b is 
given in Table 3. Assuming that GJ 1214b has a 
low albedo and its heat flow is distributed uniformly 
over its surface, its equilibriu m temperature is ~555 K 
([Miller-Ricci fc Fortnevl[2010l) . If we further assume that 
the characterization of a Rayleigh scattering tail is cor- 
rect (a = 4), then the complete data set suggests a 
mean molecular weight for GJ 1214b's atmosphere of 
around 2.7, based on the analytic model given in the Ap- 
pendix. This mean molecular weight deflnitely represents 
a hydrogen-rich atmosphere, but it leaves the problem of 
the lack of expected spectral features unresolved. 

Working from this initial prediction of a hydrogen-rich 
atmosphere. Figure [TOl shows the spectrum of a solar- 
abundance atmosphere model for GJ 1214b compared 
with the observational data. It is normalized to the ini- 
tial observation of GJ 1214b of a transit depth of 1.35% 
(1350 ppm) in the "ME arth" bandpass of 0.7-1.0 ^m 
(; Charbonneau et al.l 120091 ). This model spectrum bears 
almost no resemblance to the data, since the features are 
far too large and of different relative proportions than the 
observed features. Water features are most prominent in 
the solar model spectrum, with one that coincides with 
the suspected 1.5 /im. feature, except in magnitude. The 



Rayleigh continuum is visible as a smooth rise in transit 
depth blueward of 0.7 ^m with a similar slope to the ob- 
served rise at short wavelengths, but a far lower transit 
depth. Also, absorption dominates over scattering well 
into the optical range, so that the Rayleigh tail is pre- 
dicted primarily at wavelengths shorter than the obser- 
vations. Methane features can also be identified, such as 
those at 2.2 and 3.3 /xm, although the prominent methane 
features at shorter wavelengths overlap with water fea- 
tures. 

By ignoring one or more observed features at a time, 
the small amplitude of the remaining features can be 
explained by one or more of four possibilities: (1) the 
atmosphere contains low abundances of molecules with 
large opacities; (2) the atmosphere contains an opaque 
cloud layer that mutes the molecular absorption features; 
(3) the atmosphere contains a translucent haze layer that 
mutes the molecular absorption features; and (4) the at- 
mosphere has a high molecular weight, which results in a 
smaller scale height, a shallower atmosphere, and, thus, 
a smaller variation in transit depth. 

Figure [TT] shows the spectrum of two low-metallicity 
hydrogen-rich models. Hydrogen has a relatively low ab- 
sorption opacity compared with metal compounds, so re- 
ducing the metal abundance makes most of the features 
much smaller in amplitude. The model spectrum begins 
to conform with the data at metallicity of [Fe/H] < —2; 
the computed Rayleigh tail extends redward into the ob- 
served tail, and the spectrum is at least marginally con- 
sistent with most of the infrared observations, includ- 
ing the suggested 1.5/im feature, the i^g-band, and the 
Spitzer bands. While this is a relatively good fit, given 
the difficulty of reconciling the various observations at 
different wavelengths, the star GJ 1214 has a metallicity 
of [Fe/H] = -hO.39 (fCharbonne au et al.l[2009| ). making 
the existence of an associated planet with such a low 
metallicity extremely unlikely. We therefore do not ad- 
vance this solution. 

Figure [12] portrays spectra of model atmospheres with 
solar abundances with the addition at an opaque cloud 
layer at various pressure levels. An opaque cloud layer 
near or above the 100-mbar level ha s been suggested 
in several other models for GJ 1214b (jCroh et alj 120111: 
iBean et al.ll201ll ). and above the 10 mbar level in Berta 
et al. (2011b). We find that the spectrum becomes flat 
enough to be consistent with observations only if the 
cloud layer cuts off near or above the 1 mbar level. In 
this case, only very large features remain visible in tran- 
sit in the rarefied regions above the cloud layer. If the 
cloud tops are between 1 mbar and 0.1 mbar, the spec- 
trum is consistent with most of the measured infrared 
features (the Kg band being the primary exception) and 
would be a viable model if all of the high-amplitude fea- 
tures are rejected. Unfortunately, placing clouds at such 
a high altitude also suppresses the Rayleigh tail due to 
the molecular features, leaving the spectrum inconsis- 
tent with the observed short-wavelength rise. However, 
a translucent haze layer remains a viable alternative (see 
below). 

Figure [13] provides spectra of model atmospheres with 
high molecular weights, consisting of various proportions 
of water, with the remainder nitrogen, which we take 
to have only Rayleigh scattering opacity. With a higher 
molecular weight, the scale height of the atmosphere is 



smaller, as are the spectra features. Almost any pro- 
portion of water is consistent with most of the observed 
infrared features, again with the K^ band being the pri- 
mary exception. The best fit is an abundance of < 10%, 
water, which reproduces the suggested 1.5 /im feature in 
Berta et al. (2011b) and the K-band measurements of 
Bean et al. (2011). However, these models fail to pro- 
duce the short- wavelength rise, and the lower abundances 
fail to produce the s light observed rise in the mid-infrared 
(jDesert et al.ll201ll ). Furthermore, the molecular weights 
of water and nitrogen result in a scale height so small that 
the slope of the Rayleigh tail is much less than that of the 
observed rise. This can be confirmed using the analytic 
model in the Appendix. The slope of the Rayleigh tail 
in a pure water atmosphere (molecular weight 18) is ex- 
pected to be less than one-sixth the slope of the observed 
rise. 

As before, the addition of an opaque cloud layer only 
exacerbates the problem. This effect is illustrated in Fig- 
ure I14i specifically, with a pure water atmosphere. As 
before, the addition of clouds improves the agreement 
with the observed molecular features, but obscures the 
Rayleigh tail, leaving it a flat continuum if the cloud tops 
are above 10 mbar. 

Of the other molecules we consider, methane is the 
only one that might be consistent with the observations 
(Figure [TS]) . A pure methane atmosphere fits the large 
transit depth observed in the Kg band, but the overall fit 
is poor. A composition of < 10% methane does better, 
fitting the mid-infrared datapoints and the 1.5-/im fea- 
ture well, but it predicts additional near-infrared features 
that are not observed, and it does not fit either of the 
conflicting iiT-band and iiTg-band measurements. There- 
fore, we conclude that the principal source of absorption 
opacity is most likely water. 

4.3. Models with Haze 

On the problem of the short-wavelength rise, the high- 
molecular-weight model could possibly be salvaged with 
the addition of a haze layer. Some possible haze compo- 
nents in this temperature range are polyacetylene, tholin, 
and sulfuric acid. Sulfuric acid is known to be in the at- 
mosphere of Venus, but is considered less probable than 
carbon-based haze species, most notably due to the un- 
certainty of producing it in a water- or hydrogen-rich 
atmosphere. For this reason, and because its extinction 
opacity has a similar wavelength dependence to tholins 
at most wavelengths and particle sizes (see Figure [2]), we 
do not consider it as a solution for GJ 1214b in this work. 
However, we do produce some general models including 
sulfuric acid haze in our online library. These similari- 
ties in behavior also mean that the haze composition is 
ambiguous without more data. 

Figure [2] shows the extinction opacities of the all three 
of these haze components by particle size. At certain par- 
ticle sizes, the opacities of these particles exhibit a slope 
significantly larger than that of the Rayleigh continuum, 
meaning that they could potentially be used to fit the ob- 
served rise, even in a hydrogen-poor atmosphere. Since 
we find that hydrogen-rich models also fail to recreate 
this rise without a haze layer, we reach the conclusion 
that if the short-wavelength rise is valid, then GJ 1214b 
possesses a haze layer, which necessarily mutes molecu- 
lar absorption features, but produces a significant rise in 
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transit depth at short wavelengths. 

For a haze layer to recreate the observed rise in transit 
depth, it must be thick enough to avoid being nearly 
transparent at all wavelengths, which would produce a 
negligible effect on the spectrum. It must also be thin 
enough to avoid being nearly opaque at all wavelengths, 
since this would approach the limit of the opaque cloud 
layer, which also fails to agree with observations at short 
wavelengths. These limits vary with the particle density 
and the thickness of the haze layer. 

For polyacetylene particles of any size, the wavelength 
dependence is weak at short wavelengths. A model rep- 
resentative of the effects of monodispersed polyacetylene 
haze in a hydrogen-poor atmosphere is shown in Fig- 
ure 1161 The pressure range of the haze layer is cho- 
sen to be 0.1-0.001 mbar to maximize the contribution 
above Rayleigh scattering to total extinction at large 
radii, which should increase the slope of the scattering 
tail. This haze model has a vertical optical depth of 
r = 0.0034 at 0.85 /im, the midpoint of the normaliza- 
tion band. The small slope of the opacity curve results in 
a similarly small slope of the scattering spectrum, only 
slightly larger than the pure Rayleigh tail with no haze. 
This has a very slight effect on the rest of the spectrum. 
The resulting slope remains too small to account for the 
observed tail, so we do not support a polyacetylene haze 
in a high-molecular weight atmosphere as a model of GJ 
1214b. 

The largest slope for haze extinction opacities occurs 
with very small (< 0.03 /im) tholin particles, making this 
the best candidate to fit the observed rise in a hydrogen- 
poor atmosphere. Figure lTTl portravs several models with 
monodispersed tholin haze with the same geometric op- 
tical depth. As before, the haze layer is put between 
0.1 mbar and 0.001 mbar. At 0.85 /im, the middle of 
the normalization band, the vertical optical depth of the 
haze component in the 0.1 /im model is r = 0.068. The 
corresponding values for the 0.03 /xm and 0.01 /im haze 
models are r = 0.010 and 0.021, respectively. Changes in 
transit depth in the 0.7-1.0 /im band affect the normaliza- 
tion, but all of the models are about equally consistent 
with the infrared observations. At short wavelengths, 
the slope of the rise in transit depth is greatest with the 
smallest particles, 0.01 /im. This is consistent with obser- 
vations if the data points with the largest transit depths 
are excluded. There is some evidence for this idea in the 
qualitative bifurcation of the data blueward of 0.8 /im, 
although there are measured points from both Bean et 
al. (2011) and from de Mooij et al. (2011) at both large 
and small transit depths. 

The problem of hydrogen-rich models with hazes 
presents similar constraints. The slope required to re- 
produce the observed short-wavelength rise is a less strict 
constraint because of the lower molecular weight of the 
atmosphere, a fact that justifies a renewed interest in 
polyacetylene haze. However, a hydrogen-rich atmo- 
sphere provides an additional constraint in that large 
opacities are needed in the near- and mid-infrared to 
mute the large absorption features, especially the sug- 
gested 1.5 /im feature. 

From Figure [U the wavelength dependence of both 
haze species we consider is very similar at particle sizes 
of 1 /i,m and larger. Furthermore, the wavelength depen- 
dence is almost completely flat for 10 /xm particles and is 



relatively flat with 1 //m particles, since the particles are 
similar in size or larger than the wavelengths of light of 
interest, so their effects are roughly proportional to their 
geometric cross sections. Figure [T8l shows the effect of a 
solar- abundance atmosphere with various monodispersed 
tholin haze models with a particle size of 1 /im. At 0.85 
/im, the vertical optical depths of the models are 0.20 
for a particle density of 0.1 cm~'^ and 0.020 for a parti- 
cle density of 0.01 cm~'^. Because of the flat wavelength 
dependence, especially in the visible, the resulting spec- 
tra are flattened uniformly to varying degrees, so that 
the short-wavelength rise is suppressed, and the model 
fails to achieve the goal of fitting the short-wavelength 
data. However, with a relatively large particle density of 
0.1 cm~'^ and a high altitude of 0.1-0.001 mbar, the spec- 
trum begins to conform with the near- and mid-infrared 
data. While it is easier to keep small particles aloft, we 
do not know the degree of turbulence in GJ 1214b's at- 
mosphere, so we do not make any judgment concerning 
the plausibility of these conditions. Further research is 
needed to investigate this issue. 

Figure fTOl shows a set of models with hydrogen- rich at- 
mospheres with monodispersed tholin hazes with a parti- 
cle size of 0.1 /im. At 0.85 /im, the vertical optical depths 
of the models are 0.82 for a particle density of 1000 cm~'^ 
and 0.082 for a particle density of 100 cm~^. This figure 
demonstrates a problem that occurs in a large portion of 
the parameter space of hazes. While the haze dominates 
the atmospheric opacity at visible wavelengths, its opac- 
ity becomes too small at infrared wavelengths, resulting 
in infrared features that remain too large in amplitude to 
fit the data. This problem often persists even if the haze 
opacity saturates in the visible, and the layer becomes 
uniformly opaque at short wavelengths. Depending upon 
the choice of parameters, these models can be made to 
fit the short-wavelength rise in the data and could do so 
even if some of these data were refuted. However, in each 
case, the infrared features either are far too large to fit 
the data or are normalized to an incorrect transit depth. 

Figure [20] provides models with hydrogen-rich atmo- 
spheres with monodispersed tholin hazes with a particle 
size of 0.01 /im. At 0.85 /im, the vertical optical depths 
of the models are 0.26 for a particle density of 10^ cm~'^ 
and 0.026 for a particle density of 10^ cm~^. The results 
are much the same as for the 0.1 /im models. The opac- 
ity saturates at short wavelengths, sometimes resulting 
in a flat spectrum in the visible, while lower infrared 
opacities result in near- and mid-infrared features that 
are too large to fit the data. In other cases, the opacity 
does not saturate at these wavelengths, but the slope of 
the predicted short-wavelength tail is much larger than is 
observed. Because of the difficultly of fitting the short- 
wavelength rise and the improbability of particles this 
small, we do not advance this model as an explanation 
for the observed features. 

Polyacetylene has a spectral behavior similar to tholins 
at larger particle sizes (>1 /im), and at smaller parti- 
cle sizes (^0.1 /im), it suffers from even lower infrared 
opacities than tholins, so the same problems persist that 
occur for tholins. A haze model that might successfully 
fit both the observed short-wavelength rise in the data 
and the low amplitudes of the infrared features is one for 
which the opacity follows a Rayleigh-like A~^ behavior 
blueward of 1 /im and is roughly constant redward of 1 



fim. None of the hazes we test fits this description. 

4.4. Best-fit Models 

Our best fit models for GJ 1214b are given in Figure 
\n\ We select five models by visual inspection. While 
this is not, in principle, the most accurate method, the 
systematic problems with the data set indicate that any 
quantitative assessment of the fits would not necessarily 
be any better. Moreover, many of our theoretical spec- 
tra are similar enough that a range of parameters of the 
haze layer (in some cases as large as a factor of two) 
would produce very similar fits. Thus, any quantitative 
fit would not produce a significant gain in precision over 
our "eyeball" technique. 

Of our five selections, Model 1 is the only model we 
select that provides a good fit to the observed short- 
wavelength rise in transit depth. It utilizes a solar- 
abundance atmosphere with a tholin haze (although the 
choice of haze species is arbitrary) with a particle size of 
0.1 /im, a particle density of 100 cm~'^, and a pressure 
range of 10 to 0.1 mbar. This model has a vertical opti- 
cal depth at 0.85 fim oi t — 0.082. It is consistent with 
many of the near- and mid- infrared data points, but the 
predicted spectral features have much larger amplitudes 
than are observed. 

Models 2-4 are our best fits to the data points with 
low transit depths, with larger transit depths excluded. 
Model 2 does this with a solar-abundance atmosphere 
and a translucent tholin haze layer (the choice of haze 
species being arbitrary) with a particle size of 1 yum, a 
particle density of 0.1 cm~'^, and residing in a pressure 
range from 10-0.1 mbar; it has a vertical optical depth 
at 0.85 ^m of r = 0.20. Model 3 instead uses solar 
abundances with a uniformly opaque cloud layer, with 
cloud tops at 0.3 mbar. This model provides a better fit 
in the near- and mid-infrared than the translucent haze. 

Model 4 is our overall best fit using a high-molecular- 
weight atmosphere, specifically, 1% II2O and 99% N2 
with no clouds or haze. However, this model is not very 
sensitive to the proportion of water (within an order of 
magnitude) or to the addition of any haze that does not 
significantly obscure the absorption features. 

Model 5 is our best fit to the short-wavelength rise 
with a high-molecular- weight atmosphere. It includes 
a 1% II2O and 99% N2 atmosphere with a tholin haze 
with a particle size of 0.01 /xm, a particle density of 
lO^cm"^, and residing in a pressure range from 0.1 to 
0.001 mbar; it has a vertical optical depth at 0.85 /im of 
r = 0.021. Here, the choice of tholins is not arbitrary. 
Polyacetylene has a near-discontinuous wavelength de- 
pendence in the near-infrared at this particle size, which 
makes it unsuitable for this model. This haze model is 
relatively improbable because of its small particle size, 
and it is consistent only with the lower edge of the ob- 
served short-wavelength rise in transit depth. Therefore, 
we conclude that if this rise is valid, GJ 1214b should 
have a hydrogen-rich atmosphere. 

A hydrogen-rich model provides a better numerical fit 
to the large-transit-depth data, and the systematics pre- 
dict a low mean molecular weight for the atmosphere if 
the short-wavelength rise is valid. Furthermore, while it 
is possible to fit the low amplitudes of the observed near- 
and mid-infrared features with a hydrogen-rich model, it 



is not possible to fit the large-amplitude features with 
a hydrogen-poor model. Therefore, we conclude that a 
hydrogen-rich atmosphere is more probable overall for 
GJ 1214b, although we cannot rule out a hydrogen-poor 
model based on the current data. Because a hydrogen- 
rich model predicts larger-amplitude features, observa- 
tions of the transit depth around 2.7 and 3.3 /im, and to 
a lesser extent near 1.9 /xm, could be valuable for distin- 
guishing hydrogen-rich from hydrogen-poor models. 

We are unable to produce a model that fits all of the 
data, even after omitting either the Kg-band or K-band 
measurements. Indeed, the best fit overall is the unac- 
ceptably improbable 0.01 x solar model with no haze. If 
all of the non-conflicting data are valid, then we can say 
with confidence that there is a haze layer, but we do not 
make any predictions about its properties. Additionally, 
our methodology rules out any hazeless hydrogen-rich 
model in general and hydrogen-rich models with tholin 
hazes in the 0.01 /xm range, as well as any model in which 
the largest source of absorption is not water. We look 
forward to refining our model in light of additional data. 

5. DISCUSSION AND CONCLUSIONS 

We have developed a new capability to compute transit 
spectra of super-Earths that allows atmospheres with ar- 
bitrary proportions of common molecular species, along 
with hazes. Similar hazes have been used to characterize 
the atmospheres of giant planets, and they appear to be 
important in explaining the apparently confiicting obser- 
vations of GJ 1214b. However, additional measurements 
are needed in a range of wavelengths to determine this 
object's transit spectrum more precisely. 

If the observed short- wavelength rise is valid, then GJ 
1214b most likely possesses a hydrogen-rich atmosphere 
with a haze of small (~0.1 /im) particles, although there 
is also the possibility of an atmosphere rich in nitrogen 
and water, with very small ('^0.01 /xm), very high alti- 
tude (up to 0.001 mbar) tholin haze particles. The latter 
model is a poorer fit in the visible and is physically less 
probable, but it cannot be ruled out by the current data. 
However, the hydrogen-rich model provides a very poor 
fit redward of 1 /im. 

Alternatively, the short- wavelength measurements may 
prove inaccurate, but the near-infrared measurements 
may be accurate, including the 1.5-/im suspected water 
feature and the low-amplitude K-hsrnd measurements. 
This is a particularly plausible scenario because of the 
greater stellar variability and lower luminosity for M- 
dwarfs in the optical. In this case, in order of goodness 
of fit, GJ 1214b likely possesses (1) an N2 — H2O atmo- 
sphere, (2) a solar- abundance atmosphere, with thick, 
opaque clouds at or above the 1 mbar level, or (3) 
a solar-abundance atmosphere with a translucent haze 
with medium-to-large (>1 /xm) particles. For hydrogen- 
rich atmospheres, the models are not sensitive to the haze 
composition. High-molecular-weight models are not sen- 
sitive to the addition of hazes or clouds in general or to 
the abundance of water, within an order of magnitude. 
None of these models shows a short-wavelength rise in 
transit depth or a Ks-hand feature as deep as the one 
suggested by the data. 

In both cases, our calculations rule out a hazeless (and 
cloudless) hydrogen-rich atmosphere, which predicts a 
transit depth in the optical that is much smaller than 



the data for any reasonable metallicity. We also rule out 
a hydrogen-rich atmosphere with tholin particles in the 
'^O.Ol fiia range and a high-molecular-weight atmosphere 
with significant amounts of absorbing species other than 
water, such as NH3, CO,or CO2, which is inconsistent 
with the near-infrared data. Overall, the prospects for 
fitting the data appear better with a hydrogen-rich at- 
mosphere, but we cannot rule out a hydrogen-poor com- 
position. 

Some of the data for GJ 1214b conflict, particularly 
those in the K band and Kg band. Even for the data 
that do not directly conflict, none of the models we test 
provides a good fit to all of the observed features simulta- 
neously. However, we can rule out large classes of models 
that do not fit any of the observed features. In the case 
that all of the non-conflicting data are valid, we do not 
support any model without a haze layer. 

We propose additional observational tests to distin- 
guish hydrogen-rich from hydrogen-poor atmospheres. 
All of our hydrogen-rich models predict a strong water 
feature at 2.7 /im and a strong methane feature at 3.3 
fim, both of which should be unambiguously larger in 
amplitude than the features predicted by any reasonable 
hydrogen-poor model. The hydrogen-rich models also 
predict a water feature at 1.9 /im that could give an un- 
ambiguous result depending on the actual atmospheric 
parameters. Therefore, data at these wavelengths would 
provide strong evidence to distinguish among these mod- 
els. 

Our models correctly reproduce the expected system- 
atics of transit spectra of planets with varying metallic- 
ities, temperatures, and cloud layers. They also repro- 
duce the expected features from molecular absorption. 
Nevertheless, our models use an isothermal T-P profile. 
While this is a good first-order approximation, a more ac- 
curate T-P profile is desired to produce a more accurate 



transit spectrum, along with an atmospheric circulation 
model of the behavior at the terminator, which will not 
be in thermal or chemical equilibrium. A photochemi- 
cal model is also desired to determine the abundances of 
molecular species created or destroyed by photolysis. 

We have prepared an online library of theoretical tran- 
sit spectra produced by our calculations, sampling the 
known and expected parameter space of super-Earths. 
This library includes models with equilibrium tempera- 
tures from 300 K to 1000 K and silicate-iron planets from 
1 Afgg to 10 Mg. The library also includes models with 
atmospheric compositions of solar abundance, 0.3 x so- 
lar, 3x solar, and pure H2O, CO2, andCO atmospheres, 
as well as a variety of haze models, including sulfuric acid 
hazes. 

Very recently, we learned of the work of Murgas et al. 
(2012), who measured the transit depth of GJ 1214b at 
three additional wavelengths in the R band. For com- 
pleteness, these data points are included in Figure M 
These points are consistent with the proposed rise in 
transit depth at short wavelengths; however, their uncer- 
tainties are too large to make any definite conclusions. 
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APPENDIX 
ANALYTIC MODEL 

Lecavelier des Etangs et al. (2008) presented an analytic model of Rayleigh scattering in a planetary atmosphere as 
seen in transit. Here, we investigate the circumstances under which their model can also be applied to absorption. We 
find that the model accurately predicts the results of numerical calculations if the absorption opacity of the atmosphere 
is low, i.e., of the same order as the Rayleigh scattering opacity, but it significantly underpredicts the transit depth if 
the opacity is high. The fit is good enough at low opacities to justify this model as an estimator of the mean molecular 
weight of the atmosphere of GJ 1214b based on the slope of the short-wavelength rise in the data. 

The analytic model begins with an approximation of the optical depth of transmission at an impact parameter z, 



r(A, z) « a {\)n{z)^ 271 RpH, (Al) 

where Rp is the planetary radius, H is the scale height, and n{z) — noexp(— z/iJ) is the particle density (no being 
the density at the surface or reference altitude). Then, if i?obs is the observed radius of the planet at a particular 
wavelength, one can define r^q by -Robs = Rp + z{Tcq). That is, r^q is the optical depth at the impact parameter of the 
planet's originally observed radius. Numerical calculations suggest that for a wide range of scale heights, Teq w 0.56. 
Thus, for a given wavelength of transmitted light. 



0.56 « a(A)n(z(A))v/27ri?obs(A)i7, (A2) 

or, solving for z, the planet's effective radius lies at an altitude of. 



z(A) = i71n(Po(Tavg(A)/0.56 x ^27rRobs/kTfig), (A3) 

where CTavg is the average cross section per particle and Pq is the base pressure. 

The analytic model is valuable because it yields a single relation between the effective radius, temperature, and 
mean molecular weight of the atmosphere. If the wavelength dependence of the cross section is known, as it is for 
Rayleigh scattering, this formula allows a derivation of the temperature of the planet from its spectrum, since 

IfCTtot = (A/Ao)",then, 

U-Rohs dz 



d In A rf In A 
or, given H — kT/ fig, 



aH, (A5) 



k din A 
Alternatively, if the temperature is known, we can solve for the mean molecular weight: 

aTk din A , . „, 

M= :n^- (A7) 

g aKohs 

Plugging in the parameters for GJ 1214b yields a mean molecular weight of about 2.7. 

This model assumes that opacity is independent of temperature and pressure, which is nearly correct for Rayleigh 
scattering, but a poor approximation for molecular absorption. In the Rayleigh scattering regime, at short wavelengths 
where it dominates over absorption, we find that this model predicts the pressure at the effective radius within a few 
percent of numerical calculations. This corresponds to much less than one scale height, so the resultant transit depth 
is within a fraction of a percent of the numerical value. For absorption, the agreement is still within a factor of a 
few at most wavelengths, but the numerical calculations have a larger dispersion in effective radius. For wavelengths 
where the opacity is very high, this discrepancy can reach a factor of 100 in pressure, or about 5 scale heights. 

The dispersion and the discrepancy at the highest opacities are likely due to the analytic model's approximation 
that the opacity of the atmosphere at each wavelength is constant. Molecular absorption is often roughly constant 
over a wide pressure range, but it varies significantly at high pressures deep in the atmosphere and at the low pressures 
probed by very high opacities. Thus, the analytic model does a good job for Rayleigh scattering and when the total 
opacity is relatively low, probing higher pressures, but it becomes less accurate in cases where opacities are higher or 
are highly variable with temperature and pressure. However, it can be made accurate for particular pressure levels by 
a judicious choice of the "average" opacity. 



TABLE 1 
Properties of Transiting Earths and Super-Earths 
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Name 



Radius 

(fie) 



Mass 

(Me) 



(AU) 



Period 

(d) 



(dcg) 



GJ 1214b 


2.678 ±0.13 


6.55 ±0.98 


0.0143 


1.5803925 


88.62^0-28 


55 Cancri e^ 


2.00 ±0.14 


8.63 ±0.35 


0.01560 


0.736537 


83.4 


COROT-7b2 


1.58 ±0.1 


< 21^ 


0.0172 


0.85358 


80.1 ±0.3 


Kepler-lOb* 


1 41R+0-033 
i.4iCl_Q Q3g 


4.56ll:^^ 


0.01684 


0.837495 


84.4 


Kepler-22b5 


2.38 ±0.13 


<823 


0.849 


289.8623 


SQ 7RZ1+0-025 
89.764_(, 042 


Venus'^ 


0.950 


0.815 


0.728 


224.7 


86.61 



< 0.27 

0.17 ±0.04 





0.0068 



^ From Demory et al. (2011). 

2 From Leger et al. (2009). 

3 2cr limit. 

"' From Batalha et al. (2011). 
^ From Borucki et al. 
^ Basilevsky . 



(2011). 
Head (2003). 



TABLE 2 
Properties of Host Stars of Transiting Super-Earths-'- 



Name 


mv 


Radius 


Mass 


"^cff 


Luminosity 


[Fc/H] 


log(g) 






(Re) 


(Mq) 


(K) 


{Lq) 




(cms-2) 


GJ 1214 


14.67 


0.2110 


0.157 


3026 


0.00328 


+0.392 


4.991 


55 Cancri^ 


5.95 


0.943 


0.963 


5234 


0.582* 


+0.31 


4.45 


COROT-75 


11.668 


0.87 


0.93 


5275 


0.496 


+0.03^^ 


4.50 


Kcpler-10® 


~ 11 


1.051 


0.895 


5627 


1.004 


-0.15 


4.35 


Kcplcr-229 


11.664 


0.979 


0.970 


5518 


0.79 


-0.29 


4.44 



Uncertainties have been omitted for compactness. 

From Berta et al. (2011a). 

From Demory et al. (2011). 

From von Braun et al. (2011). 

From Leger et al. (2009). 

From Guenther et al. (2010). 

[M/H]. 

8 From Batalha et al. (2011). 

9 From Borucki et al. (2011). 
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TABLE 3 
List of Atmosphere Models for GJ 1214b 



Name 


Composition 


Haze Type 


Haze Size 
(Mm) 


Haze Density 
(cm-3) 


Haze Level F 
(mb) 


gure 


Solar 


1 X solar abundances 


None 








To 




0.1 Solar 


0.1 X solar abundances 


None 








TT 




0.01 Solar 


0.01 X solar abundances 


None 








11 




lOOmb 


Solar 


Opaque clouds 






100 


T5 




lOmb 


Solar 


Opaque clouds 






10 


~ 




Imb 


Solar 


Opaque clouds 






1 


T5 




O.lmb 


Solar 


Opaque clouds 






0.1 


12 




100h2o 


100% H2O 


None 








t3 




10h2o 


10% H2O, 90% N2 

1% H2O, 99% N2 

100% H2O 


None 

None 

Opaque clouds 








m 




Model 3 


Il3ll6ll7l21l 


100mbh2o 


100 


14 




10mbh2o 


100% H2O 


Opaque clouds 






10 


13 




lmbh2o 


100% H2O 


Opaque clouds 






1 


13 




0.1mbh2o 


100% H2O 


Opaque clouds 






0.1 


14 




100cli4 


100% CH4 


None 








15 




10ch4 


10% CH4, 90% N2 


None 








T5 




lch4 


1% CH4, 99% N2 


None 








T5 




poly 


1% H2O, 99% N2 


Polyacetylene 


0.1 


100 


0.1-0.001 


IS 




tholinO.l 


1% H2O, 99% N2 


Tholin 


0.1 


103 


0.1-0.001 


T7 




tholinO.03 


1% H2O, 99% N2 


Tholin 


0.03 


10^ 


0.1-0.001 


T7 




tholinO.Ol 


1% H2O, 99% N2 


Tholin 


0.01 


10^ 


0.1-0.001 


n 




solar 1 


Solar 


Tholin 


1.0 


0.1 


100-1 


18 




solar2 


Solar 


Tholin 


1.0 


0.1 


0.1-0.001 


18 




solar3 


Solar 


Tholin 


1.0 


0.01 


100-1 


18 




solar4 


Solar 


Tholin 


1.0 


0.01 


0.1-0.001 


Is 




solarS 


Solar 


Tholin 


0.1 


1000 


100-1 


19 




solar6 


Solar 


Tholin 


0.1 


1000 


0.1-0.001 


19 




solar? 


Solar 


Tholin 


0.1 


100 


100-1 


T9 




solarS 


Solar 


Tholin 


0.1 


100 


0.1-0.001 


T5 




solarQ 


Solar 


Tholin 


0.01 


10^ 


100-1 


20 




solarlO 


Solar 


Tholin 


0.01 


10^ 


0.1-0.001 


20 




solar 11 


Solar 
Solar 
Solar 


Tholin 
Tholin 
Tholin 


0.01 
0.01 
0.1 


10«* 
lOi^ 
100 


100-1 
0.1-0.001 [2 
10-0.1 


20 




Model 5 


0|21l 


Model 1 






Model 2 


Solar 


Tholin 


1.0 


0.1 


10-0.1 






Model 3 


Solar 


Opaque clouds 






0.3 






Model 4 


1% H2O, 99% N2 


None 








S 




Model 5 


1% H2O, 99% N2 


Tholin 


0.01 


10*5 


0.1-0.001 


El 
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Fig. 1. — Absorptive opacities of the pure molecular species H2O, CO2, H2, CH4, NH3, and CO in cm^ g~^ plotted against wavelength 
in microns at a temperature of 500 K and a total atmospheric density of 10~^ gcm~^, corresponding to a pressure of 1 mbar, which is the 
approximate pressure probed by transit spectroscopy. Vertical discontinuities extending past the bottom of the plot represent wavelengths 
at which the species in question is assumed to have zero absorption opacity. Scattering opacity is not included. 
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Fig. 2. — Extinction opactities in cm^g~^ of polyacetylene, tholin, and sulfuric acid hazes of various monodispersed particle sizes vs. 
wavelength in microns. We have assumed all of the haze species have a density of 1 gcm~^ and that they are the only contributor to 
atmospheric extinction. Thus, they are proportional to the extinction cross section of the particles multiplied by the number density and 
divided by the mass density. Constant-sloped sections indicate Rayleigh-like "^—4" behavior occurring at particle sizes significantly smaller 
than the wavelength. For polyacetylenes, this behavior occurs at long wavelengths. For sulfuric acid, it occurs at shorter wavelengths, with 
other wavelength-dependent behavior dominating in the mid-infrared. For tholins, similar spectral features dominate over Rayleigh-like 
behavior at most wavelengths. The weak oscillatory features are artifacts of the use of monodispersed particle distributions. 
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Fig. 3. — Transit depth in percent vs. wavelength in microns for generic transit spectra for hypothetical exoplanets with pure molecular 
atmospheres. The scale heights of the atmospheres are inversely proportional to their molecular weights. This proportionately affects the 
magnitudes of the spectral features, along with the slope of the Rayleigh continuum. Smooth regions of the spectra represent the Rayleigh 
continuum, with the exception of H2andC02, which additionally have smooth regions in their absorption spectra. The spectra replicate 
the molecular absorption features shown in Figure [l] Each spectrum was computed using the physical parameters and normalization for 
GJ 1214b. 
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Fig. 4. — Transit depth in percent vs. wavelength in microns for generic transit spectra computed for exoplanets with hydrogen- 
rich atmospheres of varying metaUicities. For convenience, all of the generic spectra are computed using the physical parameters (i? = 
2.678-R0, M = 6.55M0, T = 470K) of GJ 1214b, except as specified, and they are normalized to the initial mea surement of GJ 1214b. This 
normalization is 1.35% (13500 ppm) absorption averaged over the MEarth band of 0.7-1.0 fj,m (Gharbonncau et al.ll2009r) . Smooth sections 
of the graph represent the Rayleigh continuum. The continuum falls with increasing metallicity because the normalization occurs in a band 
where scattering dominates at low metaUicities, but absorption dominates at high metaUicities, with very large absorption features forcing 
the continuum down. The atmosphere is taken to be deep enough that the surface in obscured in transit at all wavelengths, and solar 
relative abundances and equilibrium chemistry are assumed. Water and methane features are prominent. Transit depth remains elevated 
between 2 and 3 ^m because of a broad collision-induced absorption feature in hydrogen. 
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Fig. 5. — Transit depth in percent vs. wavelength in microns for generic transit spectra computed for exoplanets with hydrogen-rich 
atmospheres with non-solar relative abundances. Blue: solar abundances with the abundance of methane set to zero. The resulting 
spectrum is strongly water dominated. Red: solar abundances with carbon and oxygen abundances exchanged, resulting in a C/O ratio of 
~2. This causes water features to be muted and methane features to be enhanced. The persistence of a large transit depth coinciding with 
the "water feature" at 2.7 fira is caused by a broad collision-induced absorption H2 feature. The spectra are computed using the physical 
parameters and normalization for GJ 1214b. 
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Fig. 6. — Transit depth in percent vs. wavelength in microns for generic transit spectra computed for exoplanets with solar-abundance 
atmospheres at various equilibrium temperatures. Higher temperatures result in larger scale heights and larger transit depths. The slope 
of the Rayleigh continuum also increases with temperature due to this effect. At and above 700 K, the methane abundance drops, resulting 
in weaker overlapping methane and water features (seen blueward of 1.7 ^m) and much weaker isolated methane features (seen redward 
of 1.7 /im), leaving a more water-dominated spectrum. The 4.7 fim CO line and narrow alkali metal lines begin to appear at 700 K and 
become very prominent by 1000 K. The spectra are computed using the physical parameters and normalization for GJ 1214b. 
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Fig. 7. — Transit depth in percent vs. wavelength in microns for generic transit spectra computed for exoplanets with solar-abundance 
atmospheres with achromatically opaque cloud decks at various pressure levels. The 1 bar level represents a cloudless atmosphere, that is, 
any clouds are deep enough that they are obscured by molecular absorption and scattering. Horizontal smooth regions indicate wavelengths 
where the atmosphere is transparent down to the cloud deck, so that the planet presents a fixed radius in transit with respect to wavelength. 
With very high cloud decks, many features are smoothed out, and those that remain are much weaker. If the clouds arc above the 10-mbar 
level, they also obscure the Rayleigh continuum, leaving it horizontal. The slope of the Rayleigh continuum decreases gradually before this 
point, since the clouds partially and achromatically obscure it. The spectra are computed using the physical parameters and normalization 
for GJ 1214b. 
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Fig. 8. — Transit depth in percent vs. wavelength in microns for generic transit spectra computed for exoplanets with solar-abundance 
atmospheres with tholin haze layers with different particle sizes. In each case, the haze has a number density of 100 cm""^ and extends from 
the 10 mbar pressure level to the 0.1 mbar pressure level. Horizontal smooth regions indicate wavelengths where the haze layer is opaque 
and the atmosphere is nearly transparent above it, so that the planet presents a fixed radius in transit with respect to wavelength. Sloped 
smooth regions correspond to wavelengths in which the haze layer is not opaque at the cloud tops, but it still obscures the absorption 
features. The slopes are due to the wavelength dependence of Mie scattering. The spectra are computed using the physical parameters 
and normalization for GJ 1214b. 
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Fig. 9. — Published photometric and spectroscopic transit data for GJ 1214b: transit depth in percent vs. wavelength in microns. Data 
from Sada et al. (2010) have been omitted due to large uncertainties. Data from Crossfield et al. (2011) have been omitted because they 
present a relative spectrum rather than an absolute percentage of absorption. All other published data are shown with 1-cr uncertainties. 
The horizontal "error bars" are not uncertainties, but instead represent the widths of the observational bands. 
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Fig. 10. — Transit depth in percent vs. wavelength in microns for a solar-abundance atmospheric model for GJ 1214b compared with 
the observational data. The spectrurn is normalized to the initial observation of GJ 1214b of a transit depth of 1.35% (13,500 ppni) in 
the MEarth band pass of 0.7-1.0 fim (Charbonnoau ot al. 2009). The spectral features are far too large to fit the data. Water features are 
most prominent, and one of them, at 1.5 ^m, corresponds in wavelength to one of the observed features. Distinct methane features occur 
at 2.2 and 3.3 fim, with most of the others overlapping with water features. Collision-induced absorption by hydrogen also contributes to 
the large transit depth between 2 and 3 ^m. The smooth portion of the spectrum represents the Rayleigh continuum. It has a similar slope 
to the rise at short wavelengths in the observations, but occurs only at shorter wavelengths, since it is dominated by absorption elsewhere. 
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Fig. 11. — Transit depth in percent vs. wavelength in microns for O.lx solar and O.Olx solar abundance atmosphere models for GJ 
1214b. Hydrogen has a relatively low opacity compared with metal compounds, so reducing the metal abundances results in smaller spectral 
features. At —2 < [Fe/H] < —3, the spectrum is consistent with most of the infrared observations, including the suggested 1.5 /im feature, 
the large transit depth observed in the ii's-band, and the mid-infrared Spitzer observations. Furthermore, the computed Rayleigh tail 
extends rcdward into the observed tail and has a similar transit depth to observations, unlike in the solar-abundan ce model. This is a 
relatively good fit, but we do not support it because it conflict's with GJ 1214's observed metallicity of [Fo/H] = -1-0.39 l|Charbonneau et al.l 
[2009D . 



24 



•2 1.8 

a 
s^ 
o 

CO 

^ 1.6 



CD 

o 

CD 

a. 



1.4 ^ 



1.2 



"1 1 1 1 1 1 1 1 1 T" 

Solar — abundance models for GJ 1214b with opaque clouds 



Observational data 
100 mb 
10 mb 
1 mb 
0.1 mb 




^ 0.11 



0.5 1 

Wavelength (/xm) 



Fig. 12. — Transit depth in percent vs. wavelength in microns for atmosphere models for GJ 1214b with solar abundances and an opaque 
cloud layer at various pressure levels. If the cloud tops are near or above the 1 mbar level, the spectrum is flat enough to become consistent 
with observations. In this case, only very large features should be seen in transit in the rarefied regions above the cloud layer. Specifically, 
the spectrum is consistent with most of the infrared features (the Ks band being the primary exception) . It would be a viable model for 
GJ 1214b if only low-transit-depth datapoints are considered. Unfortunately, placing clouds at such a high altitude also suppresses the 
Rayleigh tail in addition to the molecular features, leaving the spectrum inconsistent with the observed short-wavelength rise. 
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Fig. 13. — Transit depth in percent vs. wavelength in microns for atmospheric models for GJ 1214b with high molecular weights, 
specifically, various proportions of water with the remainder being nitrogen, which we take to have only Rayleigh scattering opacity. With 
a higher molecular weight, the scale height of the atmosphere is smaller, as are the amplitudes of the spectral features. Almost any 
proportion of water is consistent with most of the observed infrared features, again with the Ks band being the primary exception. The 
best fit covers a broad range of abundances of < 10%, water, which reproduces the suspected 1.5 fim feature in Berta et al. (2011b) and the 
K-hand measurements of Bean et al. (2011). However, these models fail to produce the short-wavelength rise, and the lower abundances 
fail to produce the slight observed rise in the mid-infrared. The molecular weights of water and nitrogen result in a scale height so small 
that the slope of the Rayleigh tail is much less than that of the observed short-wavelength rise. 
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Fig. 14. — Transit depth in percent vs. wavelength in microns for pure-water atmosphere models for GJ 1214b with opaque clouds at 
various pressure levels. The higher abundance of water compared to the hydrogen-rich models results in larger absorption at low pressures 
above the cloud tops and, thus, less suppression of absorption features. All of these models are consistent with the suspected 1.5 fj,m feature 
and the K-hand observations. However, the clouds suppress the Rayleigh tail, leaving it flat or nearly flat. Also, the spectra do not agree 
with the mid-infrared features. 
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Fig. 15. — Transit depth in percent vs. wavelength in microns for atmosphere models of GJ 1214b with a solar-abundance atmosphere 
with various proportions of methane with the remainder nitrogen, which we take to have only Rayleigh scattering opacity. Pure methane 
reproduces the large _K's-band observations, but also predicts features that are too large in amplitude in other near- and mid-infrared 
wavelengths. Alternatively, an abundance of <10% methane fits the suspected 1.5 fim feature well, but predicts additional near-infrared 
features that are not observed. Thus, we do not advance methane as an alternative to water to explain the observed absorption features. 



28 



1.6 



•2 1.5 

a 

o 

CO 

< 



CD 
O 

CD 



1.4 - 



1.3 - 



1.2 



T 



T 



"1 — r 



T 



T 



T 



Observational data 

No haze 

0.1 Aim, 100 cm-3 



Haze: 0.1 to 0.001 mb. 
Atmosphere: 1% HgO, 99% Ng ' 




High — molecular-weight model 

for GJ 1214b with polyacetylene haze 



_L 



_L 



J I L 



0.125 



- 0.115 



- 0.11 



0.5 1 

Wavelength (/xm) 



Fig. 16. — Transit depth in percent vs. wavelength in microns for a representative model for GJ 1214b with a 1% H2O, 99% N2 atmosphere 
with monodispersed polyacetylene haze. The pressure range of the haze layer is chosen to be 0.1-0.001 mbar to maximize the contribution 
to extinction at large radii, above the principal pressure levels probed by Rayleigh scattering, thus increasing the slope of the Rayleigh 
scattering tail. This haze model has a vertical optical depth of r = 0.0034 at 0.85 fim, the midpoint of the normalization band. The 
low slope of the Mie opacity curve (Figure O results in a similarly low slope of the resulting spectrum, only slightly larger than the pure 
Rayleigh tail without haze. There is very little effect on the rest of the spectrum. The resulting slope remains too low to account for the 
observed tail, so, if the visible-wavelength data arc vahd, we do not support a model with a polyacetylene haze in a high-molccular-wcight 
atmosphere. 
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Fig. 17. — Transit depth in percent vs. wavelength in microns for atmosphere models of GJ 1214b with a 1% H2O, 99% N2 atmosphere 
with monodispersed tholin haze. Each haze model has the same geometric depth. The haze layer is taken to be between 0.1 mbar and 
0.001 mbar. Changes in transit depth in the 0.7-1.0 ^m band affect the normalization, but all of the models are about equally consistent 
with the infrared observations. At 0.85 fim, the middle of this normalization band, the 0.1 fim haze model has a vertical optical depth of 
T = 0.068. The corresponding values for the 0.03 ^m and 0.01 fim haze models are r = 0.010 and 0.021, respectively. At short wavelengths, 
the slope of the rise is greatest with the smallest particles, 0.01 ^m, but is still less than that of the observations. This model is consistent 
with observations if the data points with the largest transit depths are excluded. 
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Fig. 18. — Transit depth in percent vs. wavelength in microns for atmospheric models of GJ 1214b with a solar-abundance atmosphere 
with monodispersed tholin haze with a particle size of 1.0 /im at various pressure levels. At 0.85 fira, the vertical optical depths of the 
models are 0.20 for a particle density of 0.1 cm~^ and 0.020 for a particle density of 0.01 cm""^. Because the wavelength dependence of the 
haze opacity is relatively flat, the amphtudes of the spectral features are uniformly suppressed. Thus, the slope of the short-wavelength 
tail is low and inconsistent with the data. For some choices of haze altitude and particle density, the near- and mid-infrared features begin 
to conform with the data. 
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Fig, 19, — Transit depth in percent vs, wavelength in microns for atmospheric models of GJ 1214b with a solar-abundance atmosphere 
with monodispersed tholin haze with particle size a 0,1 /jm at various pressure levels. At 0,85 fim, the vertical optical depths of the models 
are 0,82 for a particle density of 1000 cm~'^ and 0.082 for a particle density of 100 cm"''. These spectra can fit the observed short-wavelength 
rise in the complete data set, as well as subsets including high-depth data points, low-depth data points, and the subset of data points 
from each source, depending on the pressure level of the haze layer. However, the lower opacities in the near- and mid-infrared result in 
corresponding infrared features that are too large in amplitude to fit the data. 
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Fig. 20. — Transit depth in percent vs. wavelength in microns for atmospheric models of GJ 1214b with a solar-abundance atmosphere 
with monodispersed tholin haze with a particle size of 0.01 fim at high altitudes (between the 0.1 and 0.001 mbar levels). At 0.85 /^m, the 
vertical optical depths of the models are 0.26 for a particle density of 10 cm~ and 0.026 for a particle density of 10 cm~ . These spectra 
predict a short-wavelength rise that is steeper than in the observed data, such that the transit depth saturates near 0.5 fira. Any increase 
in particle density would increase the amount of saturation and result in a poorer fit to the data. Still, the lower opacity in the infrared 
results in infrared features that are too large in amplitude to fit the data. A decrease in haze altitude would also make the molecular 
features more prominent. Therefore, tholin particles in this size range in a hydrogen-rich atmosphere can be ruled out as an explanation 
of the spectrum of GJ 1214b. 



33 



1.8 



o 



^ 1.6 
o 

m 

< 



S 1-4 

CD 



1.2 



— I — I — I — I — I 1 r 

Best Fit Atmosphere Models for GJ 1214b 



Observational data 

Model 1 

Model 2 

Model 3 

Model 4 

Model 5 




0.13 



^eff/^s 



0.12 



zd 0.11 



0.5 1 

Wavelength (/xm) 



Fig. 21. — Transit depth in percent vs. wavelength in microns for our best fit atmosphere models of GJ 1214b. Models 1-3 use a solar- 
abundance atmosphere, while Models 4 and 5 use an atmosphere of 1% H2O and 99% N2. In Model 1, our best fit to the short-wavelength 
rise in the data, there is a tholin haze layer (although the choice of haze species is arbitrary for Models 1 and 2) with a particle size of 0.1 
jjLTCi, a particle density of 100 cm~^, and a pressure range of 10-0.1 mbar (vertical optical depth at 0.85 ^m: 0.082). Model 2 uses a tholin 
haze layer with a particle size of 1 /xm, a particle density of 0.1 cm~ , and a pressure range of 10-0.1 mbar (vertical optical depth at 0.85 
fivn: 0.20). Model 3 includes a uniformly opaque cloud layer with cloud tops at 0.3 mbar. Models 2-4 are our best fits if larger transit 
depths are excluded. Model 4 has no clouds or haze, but is not very sensitive to the addition of any haze that does not significantly obscure 
the absorption features, or to the proportion of water (within an order of magnitude). Model 5 is our best fit to the short- wavelength 
rise with a high-molecular- weight atmosphere. It includes tholin haze with a particle size of 0.01 /^m, a particle density of lO^cm"'^, and 
a pressure range of 0.1-0.001 mbar (vertical optical depth at 0.85 ^m: 0.021). However, it is consistent with only the lower edge of the 
short-wavelength data distribution and is inconsistent with many of the data points in the distribution. Therefore, we conclude that if this 
rise is valid, GJ 1214b should have a hydrogen-rich atmosphere. 



